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Abstract 

The polarization of graphene is calculated exactly within the random phase approximation for 
arbitrary frequency, wave vector, and doping. At finite doping, the static susceptibility saturates to 
a constant value for low momenta. At q = 2kp it has a discontinuity only in the second derivative. 
In the presence of a charged impurity this results in Friedel oscillations which decay with the same 
power law as the Thomas Fermi contribution, the latter being always dominant. The spin density 
oscillations in the presence of a magnetic impurity are also calculated. The dynamical polarization 
for low q and arbitrary uj is employed to calculate the dispersion relation and the decay rate of 
plasmons and acoustic phonons as a function of doping. The low screening of graphene, combined 
with the absence of a gap, leads to a significant stiffening of the longitudinal acoustic lattice 
vibrations. 

PACS numbers: 63.20.-e, 73.20.Mf, 73.21.-b 
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I. INTRODUCTION 



Recent progress in the isolation of single graphene layers has permitted the realization of 
transport and Raman experiments 1 which have stimulated an intense theoretical research 
on the properties of a monoatomic graphene sheet. Most of the unusual electronic proper- 
ties can be understood in terms of a simple tight-binding approach for the 7r-electrons of 
carbon, which yields a gapless linear band structure with a vanishing density of states at 
zero doping. 2 ' 3 

The electronic band structure of an undoped single graphene sheet allows for a descrip- 
tion of the electronic properties in terms of an effective field theory which is equivalent to 
Quantum Electrodynamics (QED) in 2+1 dimensions. Within such a framework important 
physical quantities such as the electron self-energy or the charge and spin susceptibilities 
can be calculated exploiting the special symmetries of the model. 4 

Finite doping away from half-filling qualitatively changes the above description, since 
it breaks electron-hole symmetry and also the pseudo Lorentz invariance needed for the 
equivalence with QED in two space dimensions. One is thus forced to retreat to conventional 
condensed-matter techniques such as e.g. Matsubara Green's functions. This was recently 
done for the static susceptibility of graphene at finite doping 5 . Similar calculations for bulk 
graphite had been performed by Shung et al., 6 who investigated the dielectric function and 
the plasmon behavior. Studies of static screening in graphene based on the Thomas-Fermi 
approximation had also been considered. 7 ' 8 The relation between the polarization and the 
transport properties of graphene has been recently discussed in Refs. 5,9,10. Thermo-plasma 
polaritons in graphene are discussed in citeVafek. 

In this paper, we calculate the dynamical polarization within the random phase approx- 
imation (RPA) for arbitrary wavevector, frequency, and doping. We discuss the method of 
calculation in the next section. We present in Section III results for static screening, where 
we analyze the Friedel oscillations induced by a charged or magnetic impurity, comparing 
our results with those of the two-dimensional electron gas (2DEG). Section IV discusses the 
plasmon dispersion relation and lifetime. In Section V we analyze the screening of the lon- 
gitudinal acoustic modes by the conduction electrons. Finally, Section VI presents the main 
conclusions of our work. Some of the more technical aspects of the formalism are explained 
in the Appendix. 
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II. RPA CALCULATION 



Within the effective mass approximation, and focusing on one of the two unequal K- 
points, the Hamiltonian of an hexagonal graphene sheet is given in the Bloch spinor repre- 
sentation as 2 



H = hv F Y J Hv , # k = W F Uk. (1) 




Here 4>k = k x + ik y and ipk = (ak, &k) T , where and are the destruction operators 
of the Bloch states of the two triangular sublattices. We have also introduced the Fermi 
wavevector k F which is related to the chemical potential p via kp = pjTivp. The Fermi 
velocity vp — Sat/2h is determined by the carbon-carbon distance a = 1.42 A and the 
nearest neighbor hopping energy t = 2.7 eV resulting in v F ~ 9 x 10 5 m/s. 11 We note that 
the effective Hamiltonian given above is valid only for wave vectors fc<A, where A ~ 8.25 
eV is a high-energy cutoff stemming from the discreteness of the lattice. 11 

The quantity of interest for many physical properties is the dynamical polarization, since 
it determines e.g. the effective electron-electron interaction, the Friedel oscillations and the 
plasmon and phonon spectra. In terms of the bosonic Matsubara frequencies uj n = 2im/(3, 
it is defined as 12 

P(q, iu n ) = ~ dre^(Tp(q, r)p(-q, 0)) , (2) 

where A denotes the area. The average is taken over the canonical ensemble and the density 
operator is given by the sum of the density operators of the two sub-lattices p = p a + p&. 
This amounts to working in the long- wavelength limit. To first order in the electron-electron 
interaction, we obtain 

s,s'=± 

with E ± {k) = ±hvpk — p the eigenenergies, np{E) = {er E + 1) _1 the Fermi function, 
and gs = gv — 2 the spin and valley degeneracy. A characteristic difference between the 
polarization of graphene and that of a 2DEG is the appearance of the prefactors f ss '(k, q) 
coming from the band-overlap of the wave functions 5 ' 6 

r' ( k, q) = I(i + w^±i^), ( 4) 
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where cp denotes the angle between k and q. 

At zero temperature, the Fermi functions yield simple step functions. We define the 
following retarded function by replacing iuj n — > u + i5 

9 f j2,_V^ «/ ± (k,q) 



k<D 



uj + avF^k =p |k + q|) + i5 



(5) 



where g = gs9v- The +(— ) sign corresponds to intra(inter)-band transitions and D is a 
general upper limit. 

For /i = 0, the retarded polarization thus reads 



P (1) (q,w) = -x A (q,w). 



(6) 



For /z > 0, i.e., for nonzero electron doping, the retarded polarization has an additional term 



AP (1) (q,^) =xj(q,^)+x;(q,w)- 



(7) 



In the Appendix we give additional details of the calculation of Eq. (5). Here we sum- 
marize these expressions in terms of two complex functions, F(q,uj) and G(x), defined as 

g hv 2 F q 2 



, G(x) = xVx 2 - 1 - In (x + Vx 2 - l) . 



(8) 



167T ^2 _ ^2 g 2 

From now on it is always assumed that uj > 0, noting that the polarization for u; < is 
obtained via pW(q,-w) = [P (1) (q,cj)]*. Equations (6) and (7) are then rewritten in the 
following compact form 



p( 1 »(q lW )=P (1) (q lW ) + A?( 1 )(q )W ), 



with 
and 

ApW(g,w) 



^o (1) (9,^) = -"r 



2nh 2 v 2 F 



i;2 9 ' 



-e( ^" 2 ^ + iiG 



hvFq 
hu — 2\x 



e( -2 / ,-^_ i 

hvFq 



(9) 



(10) 



,., 2/i - ftu/ 

G I — | - 27T 

nvFq 



Equations (9)-(ll) are the main result of this work. Details of the calculation are given 
in the Appendix, where we also give expressions for the real and imaginary part of the 
polarization in terms of real functions. 



Two limits of the polarization are of particular importance: (i) The long wavelength limit 
q — > with uj > vpq fixed, which is relevant for optical spectroscopy and for the plasma 
dispersion, (ii) The static case uj = with q arbitrary, which is relevant for the screening of 
charged or magnetic impurities. 

For the first case we obtain 



P^(q^0,uj) 



2 



8nhuj 



2/i 1 , 
-r— H — In 

hu 2 



2/i — huj 



2/i + huj 



- i^Q{huj - 2/i) . (12) 



And for the second case we recover previous results 5 



pm ^---^^'^^M^)' (13) 

where G<(:r) = —iG(x) is the real function (for \x\ < 1) given in Eq. (A2). Note that for 
2fi < hvpq < A the absolute value of the polarization is linear in q, and acquires rapidly 
the behavior of P^. By contrast, the polarization of the ordinary 2DEG decreases with 
increasing wave vector for q > 2k F . Furthermore, and also in contrast to the 2DEG, where 
the first derivative of the static polarization is discontinuous at q = 2k F , doped graphene 
has a continuous first derivative at hv F q = 2/i and a discontinuous second derivative. 

Table I summarizes the uj dependence of [imP^I for uj — > for a 2DEG and for doped 
graphene. While the dependence is the same for q ^ 2kp, a subtle difference appears at 
q = 2k F . This difference is intimately related to the nonanalytic behavior (as a function of 
q) of the static polarization at q = 2k F . It leads to e.g. a different power law decay of the 
electron screening, as will be shown in subsection III. Table I also indicates that, depending 
on the order in which the limits q, fi — > are taken, the low-frequency behavior may be that 
of an insulator, a metal, or a hybrid between the two. 

The polarization is a continuous function of q and uj except for the square-root divergence 
which F(q, uj) shows at uj = vpq. Figure 1 shows real and imaginary parts of the polarization 
given by Eq. (9). We note that ImP^^g,^) = for hv F q < Huj < 2fi — hvpq or < huj < 
hvpq — 2/i, and negative otherwise, while ReP^(q,uj) < for uj < vpq. 

The divergence at uj = v F q vanishes in the self-consistent RPA result of the polarization 
given by 13 

w) = i^g^y ■ <") 

Here v q = e 2 /2n q denotes the in-plane Coulomb potential in vacuum. We note that 
-Prpa(<?, vpq) = —v q so that the self-consistent polarization has a real and finite value at 



-2. ReP 1. -1. ImP 0. 




q/k F q/k F 



FIG. 1: Density plot of RePW(g,w) (left) and lmP^(q ,uj) (right) in units of n/Ti 2 v 2 F . We set 
% = vp = 1. 
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TABLE I: Low-frequency dependence of llmP^I in the limit u —> for a 2DEG and for doped 
graphene. Here y = q/2kF- 

uj = VFq- Figure 2 shows real and imaginary parts of the self-consistent polarization. While 
the singularity at u = VFq is absent, a new singularity appears in RePRPA(<?, uj) at uj oc -^/q, 
which reflects the existence of plasmons, as will be discussed in subsection IV. 



III. STATIC SCREENING 



An external charge density n ext (r) = Ze5(r) is screened by free electrons due to the 
Coulomb interaction. This results in the induced charge density 5n(r) 



8n{r) 



Ze 
4^ 



d 2 q 



Le(g,0) 



(15) 



Here e(g, 0) = lim^^o e (<L Within the RPA approximation 6 

e(q,u) = e - v q P {1) (q,uj) . 



(16) 



The effective dielectric constant eo includes high energy screening processes. We take eo 
2.4. 5 



6 




FIG. 2: Same as Fig. 1, for the renormalized polarization Prpa(q,u). 

There are two contributions to the induced charge density. A non-oscillating part comes 
from the long-wavelength behavior of the polarization. This contribution is obtained within 
the Thomas-Fermi (TF) approximation and is given by 

where a = e 2 /A-kkqUvp ~ 2.5. We note here that in a 2DEG the TF contribution also decays 
as r . 

The second contribution to the long distance behavior is oscillatory and comes from the 
non-analyticity of the polarization at Hv^q = 2/i. 14 ' 15 However, and quite importantly, in 
graphene the non-analyticity results from a discontinuity occurring only in second derivative, 
the first derivative being continuous. This leads to an oscillatory decay 

. . Zecos(2k F r) 

Sn ^ r) « M*. + 2arn ■ < 18 > 

which contrasts with the behavior of a 2DEG, where Friedel oscillations scale like 5n(r) oc 
cos(2fcpr) r~ 2 . This difference has been previously noted in Refs. 8,16, where however TF 
screening was not considered. We emphasize here that, because the TF contribution is 
of the same order of magnitude as the oscillatory part, it is essential to consider both of 
them. Furthermore we note that, while the TF contribution to screening given in Eq. (17) 
is independent of the dielectric constant 6q, the amplitude of the oscillatory part given in 
Eq. (18) decreases with increasing eo- In fact, our numerical calculations show that for 
large distances the induced density does not change sign even in the hypothetical case of 
e = where the ratio between oscillatory and TF contribution would be maximal. This 




2.5 3 3.5 

k F r / 7i 



FIG. 3: Induced charge density 5n(r) (in units of kpZe) as a function of the dimensionless variable 
kpr/n. The inset shows the r -3 decay at large distances. The high frequency screening is set 
to eo = 2.4. In the used representation the graphs are invariant under a change of the chemical 
potential. 

remarkable general property can be clearly appreciated in Fig. 3 for the particular case of 
e = 2.4: The induced density Sn(r) does not change sign, but oscillates around a finite 
offset. 

The polarization also determines the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action energy between two magnetic impurities as well as the induced spin density due 
to a magnetic impurity, both quantities being proportional to the Fourier transform of 
P^(q, 0). 17 For a magnetic impurity at r the induced spin density 8m (r) is 

1 



Sm(r) oc P (1) (r) 



4tt 2 



d 2 qP^(q,0)e i(i r . 



(19) 



Note that the magnetic response is proportional to the bare polarization, since the spin-spin 
interaction is mediated via the short ranged exchange interaction (in contrast to the long 
ranged Coulomb interaction in the case of charge response). That is why the TF contribution 
is missing and the induced spin density oscillates around zero. Figure 4 shows the numerically 
calculated pW(r) in units of k F /hvF and illustrates this behavior. Specifically in the long 
wavelength limit we obtain 



5m{r-) oc 



cos(2/c^r) 



(20) 



which is clearly seen in the inset of Fig 4. Like for the induced charge density, we find 
that the induced spin polarization 5m(r) decreases like r~ 3 for large distances. Again, this 
contrasts with the r~ 2 behavior found in a 2DEG. 17 For the particular case of undoped 
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FIG. 4: Same as Fig. 3, for the induced spin density 5m(r). 

graphene we recover the monotonous r~ 3 decay obtained in Ref. 18. 

Finally we note, that in contrast to the behavior of charge screening, where the induced 
charge scales like 1/ ji at large distances, the envelop function of the spin density is indepen- 
dent of the chemical potential at large distances. 



IV. PLASMONS 



The plasmon dispersion is determined by solving for e(q,u p — i'j) = 0, where 7 is the 
decay rate of the plasmons. For weak damping, the plasmon dispersion uj p (q) and the decay 
rate 7 are determined by 13 

Solutions to the first equation only exist for ReP^ > 0, which is the case only for finite 
doping and u > VFq. Furthermore, a stable solution requires ImPW = 0, as is the case for 
region IB of Fig. 6. Using the low-g expansion of the polarization given in Eq. (12), and 
neglecting the logarithmic correction, we obtain 

n^ M= (^iy\ (22) 

with a defined below Eq. (17). The ^fq behavior of the plasmon dispersion also appears in 
the 2DEG. 
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FIG. 5: Upper row: Solid lines show the dispersion relation for plasmons defined by Re e(q,uj) = 0. 
Dotted lines show the low-q expansion of Eq. (22), while the dashed lines represents uj = VFq. The 
crosses indicate where the plasmons acquire a finite lifetime. Lower row: Decay rate 7 of plasmons 
[see Eq. (21)] in units of chemical potential. We set h = 1. 

Outside region IB, the plasmon is damped, i.e. it has a nonzero decay rate 7. This can be 
clearly seen in Figs. 2 and 5. In the upper row of Fig. 5, we plot the exact plasmon dispersion 
for two different values of e and indicate the point at which the collective excitation becomes 
damped. The lower row shows the decay rate as obtained from Eq. (21). 

Finally, it is interesting to note that the combination of the linear dispersion relation for 
quasiparticles [electrons above (or holes below) the Fermi energy] and the plasmon dispersion 
makes it impossible for a quasiparticle of energy hu to decay into a plasmon with q < uj/vf- 
Hence, plasmons with infinite lifetime do not contribute to the lifetime of quasiparticles. 

V. ACOUSTIC PHONONS AND SOUND VELOCITY 

We now calculate the dispersion and the decay rate of acoustical phonons in graphene. 
We treat the electrons in the 7r-band of graphene as quasi-free electrons, while all the other 
electrons are assumed to be tightly bound to the carbon nuclei, thus forming effective ions 
with a positive elementary charge. In the absence of screening by the conduction electrons, 
the ions oscillate at their plasma frequency due to the long range nature of the Coulomb 
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interaction. The collective modes of the combined electron-ion plasma can be obtained 
from the zeros of the total dielectric function e to t(?, which is obtained by summing the 
contributions from ions and electrons: 19 



e tot {q, u) = e e i(<?, u) + e ion (g, u) - 1 = e - v q [P e i{q, u) + P i(m (q, w)] (23) 

Here P e i, e e i are the polarization and the dynamical dielectric function of the electrons as 
given in Eqs. (9) and (16), while Pi on ,eion are the corresponding quantities for the ions. In 
the calculation of Pi on (q,cu) we assume that the ions have a quadratic energy dispersion 
E = h 2 k 2 /2M where M denotes the ion mass. The ionic charge density is two positive 

1 /2 

charges per unit cell, so that the Fermi wave vector of the ions is k' F ~ (8n/A c ) 1 , where 
A c = 3y / 3a 2 /2 denotes the area of the hexagonal unit cell in real space. 11 This value of k' F 
is exact for field effect doping and approximate for chemical doping. 

In the calculation of P ion we assume that, for all relevant frequencies, we can take u ^> 
hk' F q/M. In the case of acoustic phonons this assumption is equivalent to v s 3> %k' F /M ~ 
1.3 x 10~Vf, where v s denotes the sound velocity. This relation is fulfilled for all meaningful 
dopings as shown below. In this regime the ion polarization is real and given by 

p ^ u) = 4^ = lKs> (24) 

where we defined E = h 2 /MA c ~7x 10 -5 eV. We note that E is of the order of the ion 
confinement energy. 

We may estimate the dispersion and the decay rate of the acoustical phonons by inserting 
Eq. (23) into Eq. (21). For finite doping /i > 0, the acoustic phonons at long wavelengths 
lie in the region 1A of Fig. 6 (defined by uj < v F q < 2fi/h — cu), where ReP^ (q, u) = 
—g/i/2irh 2 v F . In this regime the phonon dispersion is easily obtained: 

/ AnaEp \ 1/2 

uJph = —z VFq ■ (25) 

\eohvFq + gct/j, J 

To be consistent with the precondition u < vpq the expression in the square root has to be 
smaller than one. This sets a lower limit to the values of the chemical potential for which 
Eq. (25) is valid. The sound velocity v s and the decay rate 7 may be derived from Eqs. (21) 
and (25) in the limit of low q. We obtain 
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with £ = 4nE /gii. 

Some additional remarks on the validity of Eq. (25) go in place. We have already said 
that, for q — > 0, it only applies provided £ < 1. We also wish to note that the acoustical 
phonons are only well-defined if their frequency is much larger than their decay rate, i.e. if 
^ 1- Since u; p h/7 = 1 for £ = 4/5, we conclude that the notion of acoustic phonons is 
justified for £ <C 1, which corresponds to /i 3> AtiEq/ g ~ 2.2 x 1CT 4 eV. A detailed discussion 
of acoustical phonons for £ > 1 is left for future work. 

We have assumed initially that v s 3> hk' F /M ~ 1.3 x lO~ 4 v F . According to Eq. (26) this 
results in £ » 1.7 x 10~ 8 , or equivalently, /i <C 10 8 E , which is always fulfilled. 

We finish this section by estimating the sound velocity of a typical graphene sample. 
Assuming a concentration of "conduction-band" electrons of n c \ = 10 10 — 10 12 cm -2 we 
get /i ~ 1CT 2 — 1CT 1 eV. We note that this corresponds to £ ~ 2 — 20 x 1CT 3 , so that 
Eqs. (25) and (26) are applicable, resulting in v s ~ 0.05 — 0.14 -up ~ 4 — 12 x 10 4 m/s and 
7 ~ 10 -3 — 10" 2 t>i?g. These results suggest a significant enhancement of the sound velocity, 
as compared to normal metals, where v s ~ ^m/2Mv F < 10~ 2 vf- The low polarizability 
of the conduction electrons leads to a poor screening of the oscillations of the charged ions. 
On the other hand, in a semiconductor with a gap larger than the typical acoustic phonon 
frequencies, the electrons follow adiabatically the ions. In that case, the lattice vibrations 
can be described as oscillations of neutral particles. 

VI. CONCLUSIONS 

In this article we have derived a compact and closed expression for the dynamical po- 
larization of graphene within the RPA approximation. The obtained result is valid for 
arbitrary wave vector, frequency, and doping. As particular cases, we have derived the long- 
wavelength limit q — * and the static limit uj — > 0. We have employed the RPA polarization 
to calculate several physical quantities of interest in doped graphene. First we have stud- 
ied the static Friedel oscillations of the induced charge (spin) density in the presence of a 
charged(magnetic) impurity. We have found that, although the charge density does show 
oscillations around an average value, it does it without changing sign. The reason for this 
remarkable behavior is that Friedel oscillations superpose on the dominant Thomas-Fermi 
induced density, with both contributions decaying at long distances r with the same power 
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law r~ 3 . 

The dynamical polarization has been used to calculate the dispersion relation and the 
decay rate of plasmons and acoustic phonons. Like in the 2DEG case, the plasmon frequency 
shows a v /g-behavior in the long wavelength regime. We have determined the region in the 
(q, uj) plane where the plasmon is stable, as well as the decay rate in the regime where it is 
not. 

The dispersion of acoustical phonons has been shown to be strongly dependent on the 
chemical potential. In particular, we have found that the sound velocity approaches the 
Fermi velocity at low doping. However, the same limit shows an increase in the decay rate 
of acoustic phonons due to electron-hole pair excitation. 

Although we have focused on applications of the RPA calculation to the case of doped 
graphene, some aspects of the low-frequency, long-wavelength dynamics of pure graphene 
appear to be intriguing and worth studying further. 

Note added. When this work was about to be submitted, we became aware of a related 
paper. 21 Overlapping results in the two papers are in agreement. 
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APPENDIX A: CALCULATION OF THE POLARIZATION 

In the following we present some major steps of the calculation of the polarization. We 
restrict the discussion to uj > since P^\q, —uj) = [P^\q, u)] . In all the Appendix we 
set vf — h — 1, so that in the Appendix /i = fop. 
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2k F 

q 



FIG. 6: Display of the different regions characterizing the susceptibility behavior. Regions are 
limited by straight lines to = VFq (solid), to = VFq — 2fj, (dashed) and ui = 2fj, — VFq (dotted), where 
we set H = 1. 

1. Imaginary part 

The imaginary part of the functions x^(g, cu) denned in Eq. (5) has the following form 

I al3 = k / dipf(k,q)5[iU + a(k-P\k + q\)] , 
Jo 

The (^-integration yields 

' (2ak + uj) 2 - q 2 ' " 



ja/3 



q 2 — uj 2 



e(p)e( q - u>)e [k- 



q — au 



+Q(-P)e(uj - q)Q(-a) 



ei^±l-k)-e(^-k 



which is always real. The final /c-integration can now simply be performed. We obtain for 
H = 



ImP {1) (q,u) = ^-j\kJ2aI a ~(k,q,u;) 



gq 



16^u 2 -q 2 



e(u-q). (Al) 
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In order to present the result for (j, > 0, we introduce the real functions f(q, u), G>(:r), G < (x) 

J2 



9 



q 



G > (x) = xy/x 2 — 1 — cosh _1 (a;) , x > 1 , 



G<(a;) = xVl — x 2 — cos 1 (x) , \x\ < 1 



(A2) 



For the additional term at finite doping given by Eq. (7), we obtain in the language of Fig. 6 



G>(toi) -G>(^i) , 1 A 



ImAP (1) (g,cj) = — 9 - r&Va/^^g.w) = x < 

1:1 J ° t? 



7T 

-G>(^) 

-G<(^) 







1 B 

2 A 

2 B 

3 A 
3 B 



2. Real part 

The Kramers- Kronig relation valid for the retarded function KePQ 1 \q,u) reads 

>(!)/ 



ReP«(,, U ) = I/" AJ Im/ ° ^ = i£= 

W w'-w 16v/g 2 -cu 2 

— oo 



e(g-w). 



For finite doping we rewrite Eq. (7) as 



ReAP (1) (q, u) = -1L £ dkk d V J2 



2k + au + q cos if 
(k + auj) 2 — |k + q| 



(A3) 



This integral is calculated directly such that the Kramers-Kronig relation is not needed, 
here. The (^-integration yields 



ReAPd) = - 9 £ + J- 2 g jT dkJ a (k, q, u), 



(A4) 



where J a (k,q,u) is given by 

"(2aA; + cu) 2 - g 2 



J a =2tt 



w 2 — g 2 



efo- w )e(£^-*) 



+6(cu - g) 



o(n ) B(-n ) i e(^^— ^ - it) - e(fc - 
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We thus get in the language of Fig. 6 



TT 



ReAP^(q,co) = -^- + f(q,cu)x { 



-G>(^) + G>(^) 



-G<(^) + G<(te) 



■ 2/j,+u> x 



1 A 

1 B 

2 A 

2 B 

3 A 
3 B 



3. Analytic representation 

Concerning the analytic representation of the results given in Eq. (9), we note the fol- 
lowing properties (see e.g. Ref 20) 

Q(x — 1) cosh _1 (:r) = ln(x + V x 2 — 1) , 



0(1 — x 2 ) cos 1 (x) = —i ln(x + iyl — x 2 ) , 
cos _1 (— x) = 7r — cos _1 (a;) . 



Thus the functions G > (x) , G < (x) can be comprised by the single function G(x) 
x\/x 2 — 1 — ln(x + y/x 2 — 1) where 



G(x) 



G > (x) ; x > 1 

iG K (x) = -i[7T + G<(-x)] ; |x| < 1 



(A5) 



We also note the relation f(q,u) = \F(q,u)\, where f(q,u) is the real function defined 
in Eq. (A2) while F(q,u) is the complex function introduced in Eq. (8). 
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